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Asymmetric simple exclusion processes (ASEP) involve particles hopping in a 
preferred direction under hard-core repulsion that forbids double occupancy on a site 
pp. This model with periodic j21 Ej and open boundary conditions jl] as well as its 
variants involving different update schemes [Sj have been extensively studied in order 
to gain general understanding of far-from-equilibrium processes. In addition to this, 
ASEP has direct resemblance with the transport processes within the cell provided the 
dynamics of ASEP is modified by allowing attachment and detachment of particles from 
or to the environment respectively. The motor proteins which participate in cellular 
transportation by moving on linear tracks laid by long bio-molecules play the role of 
particles in ASEP El IHl E! • 

In the presence of open boundaries, one needs to think of two particle reservoirs 
attached to the boundaries which either inject or withdraw particles to or from the 
boundaries with certain specified rates. An additional reservoir is needed for the 
desorption/adsorption kinetics (Langmuir Kinetics) of particles on the lattice. The 
biased hopping of the particles, injected at one end by the reservoir, causes a finite 
current in the system even in the steady state. It is intuitively understandable that 
this particle current would help in the propagation of the boundary information to the 
bulk of the system. Thus unlike equilibrium systems, here boundaries play a crucial 
role in the steady-state dynamics and can give rise to several new features such as 
boundary-driven phase transition or production of shocks in the density profile [TIH lllj. 

The boundary- related event, that concerns us here is the appearance of localized 
shocks in the density profile. Various aspects of shocks, which are discontinuities in the 
particle density profile over a microscopic distance in the bulk, have been extensively 
studied in the past. If a and 7 are the densities maintained by the reservoirs at 
the two ends, then in the (0,7) plane there are lines a = 05(7) demarcating the 
possible phases. See Fig. H Such bulk phase diagrams are now known for many 
cases and in fact mean-field descriptions seem to give a good description of the bulk 
phase diagrams, especially for the shock formation [TIH [71 1^ ^]. Here we show the 
existence of a novel deconfinement transition of a layer near an open end as the phase 
boundary is approached from the low-density side, reminiscent of the equilibrium wetting 
transitional]. This layer with a non-bulk density profile remains attached to the end 
point, but, after deconfinement, admits the bulk density variation though with a shock. 
Let us call this special layer a shockening layer and the transition a shockening" 
transition. 

The shockening transition on the low density side is a precursor to the bulk phase 
transition, and shows power-law behaviors. This criticality is characterized by two 
length scales ^ and w, where ^ ^ 00 leads to deconfinement while w gives the length 
scale for the crossover of the surface density profile to the bulk. Though w in general 
remains finite, there is a possibility of w — > cxo which would signal a criticality at (oc, 7c) 
on the bulk phase transition line. These divergences are described by the exponents C- 
and (c defined by 

^ ~ |Aa|"^", for Aa = a — q;s(7) ^ 0— , (1) 
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Figure 1. Bulk phase boundary in the a — 7 plane, separating a low-density phase 
and a phase with a shock (insets show the density variation p vs x). The critical point 
is at (ttc, 7c)- The shock has a nonzero height on the solid line of the phase boundary, 
7 > 7c but the height vanishes on the dashed part, 7 < 7c. Arrows indicate various 
paths used in the text. 



w ~ for Aa = a — ac ^ 0—, (2) 

where Aa measures the deviation from the phase boundary for a fixed 7. In equation 
it is along a path hke path 1 in Fig. with 7 7^ 7c while for equation (j2)) it is for 
7 = 7c (path 2 in Fig. Q). By analysing a general equation for the steady state, we 
show that these two exponents C- and (c (and, in fact, several other bulk exponents 
defined below) are universal as they are determined by only a few general properties of 
the dynamics and not on details. 

A very well-studied example is the case of nonconserving ASEP of one species on a 
lattice of sites. The particles can jump to the neighboring forward site if it is empty. 
Apart from that, the bulk of the system is attached to a particle reservoir such that a 
particle can attach to (detach from) the chain with a rate uJa{uJd)- The dynamics at the 
left boundary is that of the injection of particles with a rate a while at the right it is 
withdrawal at the rate 1 — 7. Since the mean-field dynamics through average density 
like variables gives a good description of the bulk phase diagram, we adopt the same 
approach to study the shockening transition. This is expected to capture the overall 
features of the transition though fluctuations may affect the exponents. The role of 
fluctuations will be studied elsewhere. 

The mean-field equation describing the evolution of the particle density at a site i, 
is expressed in terms of the density Ui =< Ti >, where Tj is the occupation number of 
site i and (...) denotes statistical average. The dynamics is given by 

= rai_i(l - Hi) - ni{l - rij+i) + uJa{l - Ui) - LOdUi (3) 

This dynamics can obviously be extended to incorporate other effects as well. In the 
large- A^ limit, a continuum mean field approach is based on the density variable p{x) 
related to r as (Tj-ti) = p{x) 2^i^---- treating 1/A^ as the lattice spacing with 

X in the range [0, 1] . The various forms of single species dynamics studied so far can 
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be written in a general form (for the steady state) 

with fi{p), i = 0,1,2, specifying the dynamics of the system and Q = UdN. The 
steady state density profile satisfies the two boundary conditions p(x = 0) = a and 
p(x = 1) = 7 at the two ends. Here e = (2A^)~^ is the small parameter. For example, 
for the case of noninteracting particles in equation ^ the /-functions of equation (j3]) 
are 

Mp) = 1, A(p) = 2p - 1, /o(p) = K{1 -p)-p,(K = cujcu,). (5) 

The complexity of the functions /j's increases with interaction and other details but 
explicit knowledge of these functions is not essential for the analysis reported here. 

The dynamics of equation (jH) admits two special densities, (i) The Langmuir 
density corresponding to p = p^ at which /o(pl) = 0; p(x) = p^ is a particular solution 
of equation (jlj). It represents the steady-state bulk density if adsorption/desorption were 
the sole dynamics in the problem, (ii) A density p = pc at which /i(p) = 0. In case of a 
simple zero, for small deviation p = p^ + 6p, the hopping rules show a special symmetry 
of invariance of the first derivative term under 6p —6p. This is the particle-hole 
symmetry of equation (0). The bulk dynamics may not respect this symmetry in the 
general case of pc 7^ Pl- The special symmetry occurs for = Pl as e.g. for equation 
^ when K = 1. We take /2(p) 7^ for < p < 1. The cases of non-simple zeros or 
zeros of /2(p) are to be discussed elsewhere. 

To study the shockening transition, we generate a uniform approximation of the 
solution of equation Q via a leading-order boundary layer analysis ^H]- In general both 
the boundary conditions cannot be satisfied if the second derivative term (e — >■ 0) is 
ignored. As a result there appears a boundary layer at one end or a shock somewhere in 
the interior (or both). Within this special region, the second derivative term is needed. 
By neglecting appropriate terms from the original equation, one obtains two different 
solutions, to be called the outer and the inner solutions such that the outer solution is 
valid over almost the entire system and the inner solution is valid only in the region 
where the boundary layer or the shock appears. The two solutions join smoothly. As 
a general rule, if the inner solution attains a saturation then the boundary condition 
might not be satisfied by the inner solution. This is the criterion for shock formation. In 
this situation, the two boundary conditions are satisfied by two different outer solutions 
connected by the inner solution in the interior forming a shock layer. 

Without loss of generality we choose the boundary or shock layer to be at or near 
X = 1. For equation (jH) the outer and inner solutions come from 
C^Pout ^ /o(Pout) , dpi^ i^(Pm) 

dx /l(Pout) dx j2[Pin) 

where x = {x — Xd)/e is the inner variable, Xd giving the location of the layer or the 
shock, and 

np) ^ /i(po) - /i(p), with = mp). (7) 
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The matching condition Pin(5 — oo) = po = Pout(l)) for smooth joining has been 
incorporated in equation ((7j). 

Given that p(0) = a the relevant inner and outer solutions are 

Pm{x) = poSir,{x/w + and Qx = g{porit) - g{a) (left solution) (8) 

with Sin{x) — i> 1 as X — s> — oo. The functional forms of g{p) and Sin{x) depend on the 
details of /-functions. By choosing the shift = 1, is determined by 

PoSUO=l- (9) 
For example for equation (0), 

9{P) = (2p + ^ log[K - (1 + K)p]^ , (10) 

and 5*111(5;) ~ tanh(x/(2ti?) + C,), for all K. The two scales mentioned in equations (H)) 
and appear in the inner solution. There is an A^-dependence of ~ A^"''- with an 
exponent z/_ = 1 to make ^ act also as the finite size scaling variable for the transition. 
The condition for saturation of S'in(x) is 

F{p) = ioT p = ps > po. (11) 

For 7 > pin{x — > 00) the boundary condition cannot be satisfied by pm leading to 
shock formation. Therefore the phase boundary is given by 7 = Ps{po{(^s)) with shocks 
appearing for a > a<j(7). We note here that since F{p) has a vanishing derivative at 
p = pc, by Rolle's theorem of calculus, Po ^ Pc ^ Ps- In the phase with shock, the center 
of the shock , can always be chosen such that Pin(O) = pc at Xd = Xg 
Assuming simple zeros at p = po and p = p^, we write 

F{p) = -{p-p,){p-ps)<P{p), (12) 

which defines 0(p). For the the case of equation (f){p) = 1. The large x behavior is 
then given by 

dp^ P-Ps , .X . x_l /2(Ps) ..ON 

:P ~ r^, where w{a) = (p, - Po) (13) 

ax w[a) <P\Ps) 

with w depending only on a and not on 7. equation (|13|) shows w{a) as the characteristic 
length scale for approach to saturation or the bulk density, as defined earlier. A similar 
equation describes the approach to po with a scale wq (x w. This allows a practical 
definition of w as 

= - dln{p{x) - pont{x)) / dx\-^_^ , (14) 

noting that Pout(2;) is the bulk density and the deviation is only in the boundary layer 
region. 

The behavior of ^ near the phase boundary, for a fixed 7 can be determined from 
the condition at x = 0, equation 0. For 7 greater than but close to 7c, equation © 
together with equation (O, yield 

^ ~ In I7 — ps| ~ In |Aa|, or C_ = O(log), (15) 
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from equation (^Q). This scale ^ depends on both the boundary conditions. To be noted 
that the shockening transition from the low-density side determines the phase boundary. 

For a given 7, w may be made to diverge by tuning a. This locates the critical 
point on the phase boundary at (oc, 7c) such that 

ailc) - g{ac) = ^, with p^ = = p^ = 7c. (16) 

If the boundary condition at x = 1 is held fixed at 7 = 7c, then for a — > Oc, both 
Po and ps approach pc- For Po = Pc — and a = — 5a one may expand equation (jHl) 
in 5p, 5a with x = 1. Now, by definition, g'{p) = fi{p)/ fo{p), so that 

g'{p,) = 0, if pc ^ PL. (17) 

We therefore have w ^ \ps — Po|^^ ~ |Aa|~'''=, where 

Cc = 1/2, if pc 7^ PL, but Cc = 1, if Pc = PL- (18) 

The density for the critical point also shows a singular variation near the end point, 
namely, Pc — p(x) ~ y/1 — x for x — > 1 — , if Pc 7^ Pl- This follows from equation © as 
a consequence of the simple zero of /i(p). These results can be verified for the special 
case of equation (jSJ but are seen here to be of more general validity. 

The shape of the phase boundary can be determined in a similar way. For 
7 = 7c + A7, equation ISl with x = 1, on expansion gives 

-^?'(Pc)A7 + ^7"(pc)(A7)V2 + ... = g'M^a.... 

The phase boundary therefore takes the form A7 ~ lAal-^- with 

X- = l, if^?'(Pc)7^0 (e.g. i^=l) (19) 
= 1/2, if(7'(pc) = 0, (e.g. i^^l) (20) 

where K refers to equation (0). The analogous exponent x+ fo^^ 7 < 7c is discussed 
later. 

Though w measures the crossover length, right on the phase boundary, it is 
related to the height of the shock that forms beyond the phase boundary, namely, 
h = Ps — Po ^ w~^. Mean field results, based on power series expansion, seem to require 
X- < 1 so that the phase boundary is not tangential to the a-axis. This implies that 
A7 can be taken as a measure of distance r from the critical point measured along the 
phase boundary, i.e., r ~ A7. Along this curve, for 7 > 7c or r < 0, 

h ~ |r|^ with 13 = 1. (21) 

For 7 < 7c, the length w remains infinite, i.e. the shock, if formed, is of zero height. 
In fact, the boundary layer that forms does not shocken (path 3 in Fig. P). In this 
situation, there are two possibilities, either the shock height grows continuously as the 
phase boundary is crossed or the shock formation is suppressed. Though the former is 
the generic scenario, the latter situation can occur under special conditions like e.g. the 
special symmetry when pc = pi- 
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We now show that the same properties of the /-functions also determine the 
behavior in the deconfined region. The shock state, just after deconfinement, is described 
by the thickness of the deconfined layer measured by the distance of the shock from the 
boundary, A, and the height, h of the shock. If the deconfinement takes place at x = 1 
and the shock position is Xs, then A = 1 — x^. In addition the shock layer will have 
a width which we do not discuss here. The behavior of h along the phase boundary is 
also of importance, especially as one approaches the critical point. 

To discuss shock, we need the outer solution of equation ^ satisfying p(l) =7 as 

-x) = g{-i) - g{pont) (right solution) (22) 

Moreover, the shock close to the critical point is asymptotically symmetric with a scale 
w = {ps~ Po)~^ f2{Pc) /4>{Pc) as seen from equation (fT^ . For a symmetric shock, centered 
at p = pc its position Xg and height h satisfy equations (jHl) and (j221) with pout = PcTh/2 
respectively. By expanding equations (jH)), and (j221) in h and Aa = a — ac keeping 7 = 7c, 
we get (prime denoting derivatives) (path 2' in Fig. 

9\Pc)h + ^/'(Pc)/^' + ... = -g\ac)Aa + .... (23) 

A similar analysis for Xs can also be done. In the general case, pc 7^ Pl, if the third 
derivative of g is non-vanishing, equation (j23p gives 

h ~ |Aa|^', with p' = 1/3, and A = 1 - ~ |Aa|^, with C = 2/?'. (24) 

equation (j24|) defines the bulk exponents j3' and C, for h and A respectively. These 
exponents have been found for equation (jSJ with 7^ 1 [7j, but is shown here to be 
more general. There are other possibilities also. E.g. in case all the derivatives of g{p) 
vanish, no shock can exist [12], as, e.g., ioi K = 1 in equation (0). Such special cases 
will be discussed elsewhere. The exponents obtained in equation ()24|1 remain the same 
for all 7 < 7c (path 3' in Fig. ^ because in this regime the effective boundary condition 
is 7 = 7c, thanks to the formation of a nonshockening boundary layer at x = 1. A 
consequence of this is that x+ = for the shape of the phase boundary for 7 < 7c. For 
7 > 7c (path 1' in Fig. h remains 0(1) on the phase boundary so that = and 
^ = 1. It is tempting to suggest a scaling relation [3' + = 1 throughout. 

As a further example of the predictive power of this approach, we consider ASEP 
of interacting particles pBl that destroys the particle-hole symmetry. Our aim is 
to show that the same physical picture remains valid quantitatively in this case also. 
In the interior, particles at site i move to site z + 1, provided it is empty, with a rate 
that depends on the state of sites i — 1 and i + 2 as 0100— > 0010 with a rate 1 + 5 and 
1101-^ 1011 with a rate 1 — 5, while all other rates remaining the same as for equation 
(jni). For equal attachment and detachment rates {K = 1), the continuum mean-field 
equation describing the shape of the density profile is of the type equation (jH) with 
/2(p) = 1 + 5(1 — 2p), /i(p) = 1 — 2p + 5[1 — 6p(l — p)] while /o(p) remains same 
as in equation ^ with K=l, i.e. with pi = 1/2. Numerical (MATLAB) solution of 
the differential equation for relevant set of parameters actually shows that the shock is 
centered at Pc = (1 + 35 — vT+~35^)/(65), the zero of /i(p). In particular, pc < 1/2. 



Letter to the Editor 



8 



For the shockening transition, keeping 7 fixed, the relevant zero, ps, of /i(Po) — /i(p) is 
{a = 2po - 1) 

46ps = 1 + 26 -6a -{1-26^ + 26a + S^V^)^/^. (25) 

The phase boundary between the low-density and the shock phase is determined by 
the condition 7 = Ps{po, 5)- The height of the shock on the phase boundary is 7 — po 
which vanishes at the critical point 7c = pc- It follows from these that the shock near 
the critical point is asymptotically symmetric. So the critical feature, especially of the 
shockening transition of this interacting system is similar to the general case considered 
above (except for 6 = 0). This can also be explicitly verified from detail solutions. 

In summary, we have shown the existence of universal behavior associated with 
shockening transition, a precursor to the shock formation. The dynamics is described 
by a set of / functions as defined by equation (jH). The zeros pl and pc of /o(p) and 
/i(p), respectively, are the two important densities for the steady state. The shock, 
when formed, is centered around p = pc which is a point of special symmetry. The 
generic situation corresponds to Pc ^ Pl- The thickening of the layer that leads to the 
shock (called the "shockening layer") is described by two diverging scales, ^ and w whose 
critical behaviors are determined by the nature of fi{p) around p = pc. Furthermore, our 
analysis and results suggest the possibility of more complex dynamics involving different 
symmetries, though such model systems are not known till now. 

SM acknowledges financial support from IIT, Kanpur and thanks Bhaskar Dasgupta 
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